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Abstract. Two discretizations, linear and nonlinear, of basic notions of the complex analysis are 
considered. The underlying lattice is an arbitrary quasicrystallic rhombic tiling of a plane. The linear 
theory is based on the discrete Cauchy-Riemann equations, the nonlinear one is based on the notion of 
circle patterns. We clarify the role of the rhombic condition in both theories: under this condition the 
corresponding equations are integrable (in the sense of 3D consistency, which yields also the existense of zero 
curvature representations, Backlund transformations etc.). We demonstrate that in some precise sense the 
linear theory is a linearization of the nonlinear one: the tangent space to a set of integrable circle patterns 
at an isoradial point consists of discrete holomorphic functions which take real (imaginary) values on two 
sublattices. We extend solutions of the basic equations of both theories to Z"*, where d is the number of 
different edge slopes of the quasicrystallic tiling. In the linear theory, we give an integral representation 
of an arbitrary discrete holomorphic function, thus proving the density of discrete exponential functions. 
We introduce the d-dimensional discrete logarithmic function which is a generalization of Kenyon's discrete 
Green's function, and uncover several new properties of this function. We prove that it is an isomonodromic 
solution of the discrete Cauchy-Riemann equations, and that it is a tangent vector to the space of integrable 
circle patterns along the family of isomonodromic discrete power functions. 

1 Introduction 

There is currently much interest in finding discrete counterparts of various structures of the classical 
(continuous, smooth) mathematics. In the present paper we are dealing with the discretization of 
the classical complex analysis. 

There are two approaches to this problem. The first one, which we shall call the linear theory, 
is based on a discretization of the Cauchy-Riemann equations. Since the latter are linear, straight- 
forward discretizations are linear as well. A discretization preserving apparently the most number 
of important structural features has been developed in |FI IDll lD2l IMll IK] . The first two references 
are dealing with discrete holomorphic functions / : — > C on the regular square lattice, satisfying 
the following discrete Cauchy-Riemann equations: 
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A pioneering step was undertaken by Duffin |D2j . where the combinatorics of 1? was given up in 
favor of arbitrary planar graphs with rhombic faces. A far reaching generahzation of these ideas is 
given in |Mlj . where the hnear theory is extended to discrete Riemann surfaces. Planar graphs with 
rhombic faces are called critical in |Mlj . Kenyon developed a theory of the Dirac operator and 
constructed Green's function in the framework of the linear theory on critical graphs. See |CY1 lU] 
for combinatorial, resp. numerical aspects of Green's functions on graphs. 

The second approach, which we will call the nonlinear theory^ is based on the ideas by Thurston 
[T] . and declares circle patterns to be natural discrete analogs of analytic functions jBeStlDSllSchl lS]. 
One of the most important achievements of this theory is the proof that the Riemann map can 
be (constructively) approximated by circle packings |RS[ IMRl IHSj . The variational approach to 
circle patterns is discussed in detail in |BSp| . The word "nonlinear" refers to the basic feature of 
equations describing circle patterns. Often, the so-called cross-ratio system is used for this. For a 
function / : ^ C on the regular square lattice, this system was introduced in [NT 
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(/m+l,ra fm,n){fm+l,n+l fm,n+l) 

(/m,n+l fm,n){fm+l,n+l fm+l,n) 

For circle patterns with more sophisticated combinatorics, a generalization of this system to an 
arbitrary quad-graph (planar graph with quadrilateral faces) is required |BSj . 

It is not difficult to see in what sense solutions of equations like Q can be considered as 
discretized analytic functions. Indeed, assume that is embedded in the complex plane C with 
the grid size e, i.e., the pair (m, n) G Z"^ corresponds to {m + in)e E C. Then restrictions of analytic 
functions to this grid satisfy the corresponding equations up to O(e^). More precisely, if / : C ^ C 
is analytic, then 

f{z + s + ie) - f{z) 

and 

{f{z + e)- f{z)) {f{z + e + ie)- f{z + ie)) 
[fiz + ie) - f{z)) {fiz + e + ie)- f{z + e)) 

Similar relations hold on more general graphs. 

For a long time, the linear and the nonlinear theories of discrete analytic functions were consid- 
ered separately. In the present paper, we show that in some precise sense the former is a linearization 
of the latter. We work in the set-up of rhombic tilings of a plane. The theory becomes especially 
rich for quasicrystallic tilings, - those with a finite number of different edge slopes. This class 
includes double periodic tilings (which are naturally considered on a torus), as well as non-periodic 
ones, like the Penrose tiling. We clarify the importance of rhombic embeddings of quad-graphs in 
both the linear and the nonlinear theories. Namely, we show that the rhombic property implies 
(actually, is almost synonymous with) integrahility. Note that interrelations of circle patterns with 
the theory of integrable systems were already uncovered and studied in |BP| lABU IAB2| IBHS| IBH) . 
Note also that some of the ideas behind our unified treatment of integrability of linear and non- 
linear systems, such as the use of zero curvature representations in both situations, are similar to 
the philosophy of Fokas's unified transform method for linear and nonlinear differential equations 
based on the Riemann-Hilbert boundary problem |Foj . Our main results are the following. 



Discrete Cauchy-Riemann equations on a rhombically embedded quad-graph P, with weights 
given by quotients of diagonals of the corresponding rhombi, are integrable. Integrability is 
understood here as 3D consistency |BSj . Therefore, discrete holomorphic functions on rhombic 
embeddings can (and should) be extended to multidimensional lattices. In particular, discrete 
holomorphic functions on a quasicrystallic rhombic embedding T) with d different edge slopes 
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can be considered as restrictions of discrete holomorphic functions on Z to a certain two- 
dimensional subcomplex Qx> in 



• Cross-ratio equations on a rhombically embedded quad-graph V, with cross-ratios read off the 
corresponding rhombi, are integrable as well. Therefore, solutions of the cross-ratio equations 
on a quasicrystallic rhombic embedding V are naturally extended to Z'^. 

• For a circle pattern, the centers and the intersection points of the circles yield a solution 
of cross-ratio equations, with the cross-ratios depending on the pairwise intersection angles 
of the circles. We say that a circle pattern is integrable, if the corresponding cross-ratio 
system is integrable. The combinatorics and intersection angles belong to an integrable circle 
pattern, if and only if they admit an isoradial realization. This latter realization gives a 
rhombic immersion of the corresponding quad-graph, and generates also a dual isoradial 
circle pattern. An integrable circle pattern can be alternatively described by the radii of the 
circles and the rotation angles of the configurations at the intersection points with respect to 
the isoradial realization. These data comprise a solution of an integrable Hirota system. 

• The tangent space to the set of integrable circle patterns, at the point corresponding to an 
isoradial pattern, coincides with the space of discrete holomorphic functions on the corre- 
sponding rhombically embedded quad-graph, which take real (resp. pure imaginary) values 
on the white (resp. black) vertices. This holds in the description of circle patterns in terms 
of circle radii and rotation angles at the intersection points (Hirota equations). Discrete 
holomorphic functions obtained from these ones by discrete integration, comprise the tan- 
gent space to the set of integrable circle patterns, described in terms of circle centers and 
intersection points (cross-ratio equations). 

• We define (in the linear theory) discrete exponential functions on Z'^, and prove that they are 
dense in the space of discrete holomorphic functions, growing not faster than exponentially. 

• We define (in the linear theory) a discrete logarithmic function on Z'^, or, better, on a branched 
covering of certain d-dimensional octants ^ Sm C Z*^, m = 1, . . . , 2d. On each such octant, the 
discrete logarithmic function is discrete holomorphic, with the distinctive property of being 
isomonodromic, in the sense of the integrable systems theory. We show that the real part of 
the discrete logarithmic function restricted to a surface Jlx> in coming from a quasicrystallic 
quad-graph T> is nothing but Green's function found in pj. The integral representation of 
Green's function given in [K] is derived within the isomonodromic approach. 

• We define (in the nonlinear theory) discrete power functions w^^^ (resp. z'^) on the same 
branched covering of octants Sm <Z Z'^ , m = 1, . . . , 2d, where the discrete logarithmic function 
is defined. On each such sector, discrete w^~^ (resp. z'^) is an isomonodromic solution of 
the Hirota (resp. cross-ratio) system. The tangent vector to the space of integrable circle 
patterns along the curve consisting of patterns w'^~'^ , at the isoradial point corresponding to 
7 = 1, is shown to be the discrete logarithmic function. 

In conclusion, we point out some generalizations of the concepts and results of this paper for the 
non-rhombic case. 

Acknowledgements. Numerous discussions and collaboration with Boris Springborn were 
very important for this research. We thank also Tim Hoffmann, Ulrich Pinkall and Giinter Ziegler 
for discussions. 

^We use this term for a subset of Z"^ defined by fixing one of 2'' possible combinations of signs of the coordinates. 
An octant in the proper sense corresponds to d = 3, while by d = 2 this object is called quadrant. 
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2 Discrete harmonic and discrete holomorphic functions on graphs 



We denote by V{G), E{Q) and E{Q) the sets of vertices, undirected and directed edges of a graph 
Q, respectively. Let there be given a complex-valued function : E{Q) ^ C on the edges. Then the 
Laplacian A corresponding to the weight function v is the operator acting on functions / : V{G) C 
by 

(A/)(xo) = ^(^0, x)if{x) - /(xo)). (3) 

Here the summation is extended over the set of all vertices x connected to xq by an edge. We will 
use the notation star(2;o) = star(xo; G) for the set of all edges of Q incident to xq, see Fig.^ 

Definition 1 A function f : V{G) ^ C is called discrete harmonic (with respect to the weights 
v), ifAf = 0. 

Of course, the most interesting case of these notions is that of real positive weights i' : E{Q) — s- K+. 



X2 




Figure 1: The star of the vertex xq in the graph Q . 

In the continuous case, there is a canonical correspondence between harmonic and holomorphic 
functions on C: the real and the imaginary parts of a holomorphic function are harmonic, and any 
real-valued harmonic function can be considered as a real part of a holomorphic function. This 
relation can be generalized for functions on graphs, but these two classes of functions live then on 
different graphs. Discrete holomorphic functions live on quad-graphs. 

Definition 2 A cell decomposition T> of the plane C is called a quad-graph, if all its faces are 
quadrilaterals. 

A more general version of this definition deals with cell decompositions of an arbitrary oriented 
surface. So, quad-graphs are not just graphs, but are additionally assumed to be embedded in an 
oriented surface; we will deal with the case of C only. 

To establish a relation with discrete harmonic functions, we consider the latter ones on graphs 
Q with an additional structure, namely on those that come from general (not necessarily quadri- 
lateral) cell decompositions of C. We will denote by F{Q) the set of faces (2-cells) of Q. To any 
such Q there corresponds canonically a combinatorial quad-graph called its double (or diamond 
|Mlj '). constructed from Q and its dual Q* . Recall that, in general, a dual cell decomposition Q* 
is only defined up to isotopy, but it can be fixed uniquely with the help of the Voronoi/Delaunay 
construction. The dual Q* is characterized as follows. Vertices of Q* are in a one-to-one cor- 
respondence to faces of Q, see Fig. |2I Each e E E[Q) separates two faces of which in turn 
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correspond to two vertices of Q* . It is declared that these two vertices are connected by the edge 
e* G E{Q*) dual to e. Finally, the faces of Q* are in a one-to-one correspondence with the ver- 
tices of Q: if xo € y{G)-: and xi, . . . ,x„ G V{Q) are its neighbors connected with xq by the edges 
^1 = {xq,xi), . . . , e„ = (xo,x„) G E{Q), then the face of Q* corresponding to xq is defined by its 
boundary ej U . . . U e* (cf. Fig. If one assigns a direction to an edge e G E{Q), then it will be 
assumed that the dual edge e* G E{Q*) is also directed, in a way consistent with the orientation 
of the underlying surface, namely so that the pair (e, e*) is oriented directly at its crossing point. 
This orientation convention implies that e** = — e. 




Xti 

Figure 2: Vertex of Q* dual to a face of Q. Figure 3: Face of Q* dual to a vertex of Q. 

Now the double T> is constructed from Q* as follows. The set of vertices of the double T) 
is V{V) = V{g) U V{g*). Each pair of dual edges, say e = (xo,xi) G E{g) and e* = (yo,yi) e 
E{Q*), defines a quadrilateral {xo,yo,xi,yi). These quadrilaterals constitute the faces of the cell 
decomposition (quad-graph) V, see Fig.^ The edges of V belong neither to E{g) nor to E{Q*). 
A star of a vertex xq G V{g) produces a flower of adjacent quadrilaterals from F[T)) around the 
common vertex xq, see Fig.|21 

Observe that the double T> is automatically bipartite, since its vertices V{T>) are decomposed 
into two complementary halves, ViT^) = V{Q) U V{Q*) ("black" and "white" vertices), such that 
the ends of each edge from E{T)) are of different colours. An arbitrary quad-graph embedded 
in C is automatically bipartite, and the above construction can be reversed for it, to produce a 
cell decomposition Q along with its dual Q* . The decomposition of V{T>) into V{Q) and V{Q*) is 
unique, up to interchanging the roles of Q and Q* . Edges of Q (say) connect two "black" vertices 
along the diagonal of each face of T>. 




Figure 4: A face of the double T) Figure 5: Faces of T> around vertex xq. 
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Let there be given a function v : E{Q) ^ C on undirected edges of Q. (It is assumed that 
both directed representatives ite of any edge carry the same value = e) as the underlying 
undirected one.) Extend the function u to undirected edges of Q* according to the rule 

= 1/Ke). (4) 



Definition 3 A function f : V{'D) — > C is called discrete holomorphic (with respect to the 
weights v), if for any positively oriented quadrilateral {xo,yo,xi,yi) G F{I)) there holds: 

—— — — = ii^{xo, xi) = — -. (5) 

/(xi)-/(xo) iv(yo,yi) 

These equations are called discrete Cauchy-Riemann equations. 

Again, the most interesting case corresponds to the real positive weights v : E{Q) U E{Q*) R+. 
The theory of discrete holomorphic functions was developed in jMlj . In ^ a discrete Dirac operator 
was introduced, the kernel of which consists of discrete holomorphic functions. In the present paper, 
like in |Mlj . our attention belongs not to the discrete Dirac operator but to its kernel only. The 
next statement follows immediately. 

Proposition 4 a) If a function f : V{T)) ^ C is discrete holomorphic, then its restrictions to 
V{Q) and to V{G*) are discrete harmonic. 

b) Conversely, any discrete harmonic function f : V{Q) — > C admits a family of discrete 
holomorphic extensions to V{V), differing by an additive constant on V{Q*). Such an extension is 
uniquely defined by a value at one arbitrary vertex y £ V{Q*). 



3 Rhombic embeddings and labelings of a quad-graph 

The paper |KSj studies rhombic embeddings of a quad-graph V in C, i.e., embeddings with the 
property that each face of V is mapped to a rhombus. A combinatorial criterion for the existence 
of a rhombic embedding of a given quad-graph V found in |KSj is as follows. 

Definition 5 A strip S in D is a path {aj}j%_^ in D* with the following property: for any two 
consecutive dual edges a*, a*_^_l € S C E{V*) with the common point a*^a*_^_-^ = qj G ¥{!)*) ~ F{'D), 
the corresponding edges aj,aj+i € E{T>) are two opposite sides of the quadrilateral qj. The edges 
{aj}'jl_^ are called the traverse edges of the strip S. 

Theorem 6 |KSj A planar quad-graph T> admits a rhombic embedding in C if and only if the 
following two conditions are satisfied: 

• No strip crosses itself or is periodic. 

• Two distinct strips cross each other at most once. 

A rhombic embedding determines rhombus angles that are naturally assigned to the edges of G 
and G* , see Fig.lHI Such systems of rhombus angles (j) : E{Q) U E{Q*) (0, tt) are characterized, 
according to IKSj . by the following two conditions: first, 

cl){t*)=7T-ct>{c), yceEiG), (6) 
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Figure 6: A rhombic embedding of a quadrilateral {xo,yo,xi, zi) € F(V), e = {xo,xi), e* = {yo,yi) 



and, second, 

</>(e) = 27r, J2 '^(^*) = ^ ^(^)' ^ ^(^*)- 

eestar(xo;5) c*estar(j/o;C?*) 

As mentioned in |KSj , to each rhombic embedding of T> there corresponds a set of parallelogram 
embeddings (wherein each face is mapped to a parallelogram), which are obtained by replacing each 
traverse edge of a strip with a real multiple (a different multiple for each strip). 

Definition 7 A labeling is a function a : E{V) — > C such that a{—a) = — Q(a) for any a G E{'D), 
and the values of a on two opposite and equally directed edges of any quadrilateral from F{'D) are 
equal to one another. 




Figure 7: Labeling of directed edges Figure 8: Labeling of undirected edges 

This definition is illustrated in Fig.d Note that if edges of any given face of T) are directed 
as on this figure (from black to white), then any two opposite edges carry opposite labels. For 
any labeling a : E{T>) — > C of directed edges, the function can be considered as a labeling of 
undirected edges, i.e., as a function : EiV) —>■ C such that its values on two opposite edges of any 
quadrilateral from F{T>) are equal to one another, see Fig.|HI Conversely, any labeling of undirected 
edges comes as a square of some labeling of directed edges. 

The existence of a labeling a : E('D) — > C is equivalent to the existence of a parallelogram 
immersion of the quad-graph D. Indeed, given a parallelogram immersion p : V(T>) C, one 
defines canonically a labeling by 

a{x, y) = p{y) - p{x), V(x, y) G E{V). (8) 

Conversely, given a labeling a : E{T)) — > C, the formula correctly defines a function p : V{T>) — > 
C and assures that the p-image of any quadrilateral face of P is a parallelogram. If the labels a 
take values in = {0 G C : \9\ = 1}, then the corresponding immersion is rhombic. 
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Definition 8 A parallelogram immersion p : V(T>) — > C o/ a quad-graph V is called quasicrys- 
tallic, if the set of values of the corresponding labeling a : E{'D) C, defined by is finite, say 
A = {±ai, . . . ,±arf}. 

It will be supposed that any two non-opposite elements of A are linearly independent over M. This 
implies, in particular, that all parallelograms are non-degenerate. 

It will be of a central importance for us that any quasicrystallic parallelogram immersion p can 
be seen as a sort of a projection of a certain two-dimensional subcomplex (combinatorial surface) 

of a multi-dimensional regular square lattice Z*^. The vertices of Q.ti are given by a map 
P : V{T)) — i- constructed as follows. Fix some xq € V{T>)^ and set P{xq) = 0. For all other 
vertices of their images in are defined recurrently by the property: 

• For any two neighbors x, y G if p{y) — p{x) = itaj E A, then P{y) — P{x) = ztej , 

where is the i-ih coordinate vector of Z'^. The edges and faces of f]© correspond to edges and 
faces of T). So, the combinatorics of is that of P, and therefore Theorem El can be used to 
decide whether a given two-dimensional subcomplex of Z'^ corresponds in this way to some rhombic 
embedding of a quad-graph in C. 

4 3D consistency 

We now study a question about integrability of the discrete Cauchy-Riemann equations ©. These 
equations are just a specific linear issue of general equations on quad- graphs jBSj 

^(/(a:o),/(yo),/(xi),/(yi)) =0, (9) 

relating four fields / sitting on the four vertices of an arbitrary (oriented) face (xq, yo; xi,yi) G F{T>) 
of a quad-graph T>. Here the function $ may depend on some parameters (in the case of discrete 
Cauchy-Riemann equations these are the weights v), and it is supposed that equation ® is uniquely 
solvable for any one of the fields in terms of other three (which is, of course, the case for discrete 
Cauchy-Riemann equations with non- vanishing weights v). 

The approach pushed forward in [BS j is based on the idea that integrability of such equations 
on quad-graphs is synonymous with their 3D consistency. To describe the latter notion, we extend 
the planar quad-graph T> into the third dimension. Formally speaking, we consider the second 
copy T> oiT> and add edges connecting each vertex x G Vi'D) with its copy x € V{'D). On this way 
we obtain a "three-dimensional quad-graph" D, whose set of vertices is 

y(D) = V{V)UV{V), 

whose set of edges is 

-E(D) = E{V) U E{V) U {{x, J) : X G V{V)}, 

and whose set of faces is 

F(D) = F{V) U F{V) U {(x, y, y, x) : (x, y) G E{V)]. 

Elementary building blocks of D are cubes (xq, yQ,xi,yi,XQ,yQ,xi,yi), as shown on Fig.l^l 

Clearly, if T) is bipartite, then so is D: each x G V{T>) has the colour opposite to the colour of 
its counterpart x G V{'D). Hence, we can extend the "black" graph ^ to a 3D object G, with edges 

E{G) = E{g) U E{g^*) U {(x, y) : X G V{g), y G V{g*), (x, y) G E{V)}. 

Edges of E{G) within Fig.lHlare (xo,xi), (yo,yi), {xo,m), (xo,yi), (xi,yo), and (xi,yi), forming 
the black tetrahedron. Similarly, we have a 3D white graph G*. 
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Figure 9: Elementary cube of D 



Definition 9 Equation ^ is called 3D consistent if it can be imposed on all faces of any ele- 
mentary cube ofD, in such a manner that opposite faces carry one and the same equation (i.e., 
the same parameters). 

This should be understood as follows. Consider an elementary cube of D, as on Fig.|21 Suppose 
that the values of the function / are given at the vertex xq and at its three neighbors yo, yi, and 
xq. Then equation Q uniquely determines the values of / at xi, yo, and yi. After that equation 
© delivers three a priori different values for the value of the field / at the vertex xi , coming from 
the faces {yo, xi,xi, yo), {xi,yi, yi,xi), and (xq, yo, xi, yi), respectively. The 3D consistency means 
that these three values for f{xi) actually coincide, independently on the choice of initial conditions. 

As discussed in detail in |BSj , the 3D consistency of a given system allows one to construct 
Backlund transformations and to find in an algorithmic way a zero curvature representation for it, 
which are traditionally considered as main attributes of integrability. Briefly, the constructions are 
as follows. 

1) Given a solution / : ^ C to © and an arbitrary value f{xo) = fo at some vertex 
Xq € V, the 3D consistency allows one to extend the solution / to the whole of V^D). Its restriction 
to V(T>) is thus a well-defined function / : 1^(1)) C which also solves the original equation ©. 
Setting /(x) = f(x) for all x € V{'D), one can interpret this function as / : V(T>) — > C, and this / 
is called the Backlund transformation of / (defined by the value fo and the parameters sitting 
on the vertical faces). 

2) Suppose that the function <I>(ui, n2, M3, M4) in Q is affine-linear in all its arguments, so 
that this equation can be solved uniquely for an arbitrary argument Ui in terms of other three 
arguments, the solution being given by a fractional-linear function. For an arbitrary edge a = 
(x, y) £ E{T>), consider the vertical face (x, y, y, x) € F{0) over this edge. The solution of equation 
^{f{x), f{y), f{y), f{x)) = can be written as 

f{y) = L{f{y),f{x))- fix), (10) 

where L{f{x),f(y)) £ PGL2{C), and the standard notation for the action of PGL2{C) on C by 
Mobius transformations is used: 

/a b\ _ au + b 
\c dj cu + d' 

One assigns the matrix above to the edge a, so that L{a) = L{f[y), f{x)). Now it follows from the 
3D consistency that for an arbitrary face (xq, yoi ^^i, yi) € F{'D) one has: 

L{f{xi), /(yo)) L{f{yo),f[xo)) = L{f{xi),f{yi)) L(/(yi), /(xq)). (11) 

This expresses the flatness of the discrete connection L on D with values in PGL2{C), hence 

is called the zero curvature representation of system @. It is often possible to use suitable 

normalizations in order to lift this representation to the one with values in GL2(C). 
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5 3D consistent Cauchy-Riemann equations 



To apply the notion of the 3D consistency to the discrete Cauchy-Riemann equation one has 
to explain how to impose it on the vertical faces of D. For this, we assume that the function u is 
extended to E(G) U E{G*), still satisfying the condition z^(e*) = l/z^(e), and with an additional 
condition that opposite edges carry the same values of v. 

An interesting problem is, of course, to find functions v on the "ground floor" E{Q) U E{Q*) 
which can be extended to the edges of E{G) U E{G*) lying in the "vertical" faces to give a 3D 
consistent system. 

Theorem 10 The function v : E{g) U E{g*) C can be extended to E{G) U E{G*) giving a 3D 
consistent system of discrete Cauchy-Riemann equations, if and only if the following condition is 
satisfied: 

n 1^=1. n ^^=1. v.„ . » e v(e-). (12) 

eestar(a;o;6) c*estar(yo;5*) 

Proof. Consider a flower of quadrilaterals around xq, with = {xo,Xk), = {yk-i^Uk) (in 
notations of Fig.[^. Build the extension of this flower to the third dimension (for one of its petals 
corresponding to A; = 1 this extension is shown on Fig.|^. Denote 

v{yk-i,yk) = = Vk, v{yk,Xo) = l^k- (13) 



Lemma 11 Discrete Cauchy-Riemann equations are 3D consistent on the cube over the k-th petal, 
if and only if 

1 + Vk^ik-i - T^kl^k + fJ'k-ilJ'k = 0. (14) 

Proof of Lemma II 11 Consider the elementary cube on Fig.|^ corresponding to A; = 1. On the 
first step of checking the 3D consistency we find: 

/(xi) = f{xo)+iui{f{yo)-f{yi)), 
fim) = f{xo) + iMfivo) - f{xo)), 
fim) = /(xo)+i/ii(/(yi)-/(xo)). 

On the second step we find (from the condition that opposite faces support the same equations): 

/(xi) = f{xo) + iui{f{yo)-f{yi)), 

= fiyi) + iMfixi) - fim)), 
= /(yo) + w(/(2;i) - /(yo))- 

After simple computations we find: 

/(xi) = {l + i^ifJ-o-i^iHi)f{xo)-uifiof{yo)-iriyiHif{yi), 

= -/^ow/(2o) - z^i/^o/(yo) + (1 + 1^1/^0 + /^om)/(yi), 
= -/^om/(^o) + (1 - J^im + ^om)/(yo) + ^ii^ifiyi)- 

Comparison of these expressions leads to 1 + uifiQ — vi^i + /^o^i = 0, which proves the lemma. I 
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Continuing the proof of Theorem IIUI we derive from 1)141) : 

where the standard notation for the action of PGL2{C) on C by Mobius transformations is used. 
Starting with an arbitrary /iQ, we can define all /i/^'s consecutively. This procedure is consistent, if 
running around xq returns the value of fiQ we started with. This holds for any /iQ, if and only if 

the matrix product (^^\ ^ scalar matrix. It is easy to see by induction that the above 

matrix product may be presented as 



A B 
-B A 



^ with A = ]^{j\{vk + i)+J[{T^k-i)), B = j.{j\{vk + i)-\{{vk-i) 



k k k 

Therefore, a necessary and sufficient condition for this matrix to be scalar is 

I'k + i 



which is equivalent to the first equality in ^Vl\ . because of Uk = ^{'^k) ~ The second 

condition in ()12j) is proved similarly, by considering a flower of quadrilaterals around yo £ V{G*). 

m 

As pointed out above, the most interesting case is when u takes values in K^-. In this case we 
will use the notation 

i/(e)=tan^, 0(e) g(0,^). (15) 

The condition z^(e*) = l/j^(e) is translated in this case into The integrability condition H12|) 
takes in this case the form 

Yl exp(2<A(e)) = 1, J] exp(i0(e*)) = 1, Vxq £ V{g), yo € V{g*)- (16) 

egstar(a;o;(7) c*gstar(j/o;G*) 

The latter condition is a generalization of ((TJ, and is equivalent to saying that the system of angles 
(j) : E{Q) U E{Q*) — > (0, vr) comes from a realization of the quad-graph P by a rhombic ramified 
embedding in C. Flowers of such an embedding can wind around its vertices more than once. 

Lemma 12 Let a quad-graph V be a double for a pair of dual cell decompositions Q, Q* . Let 
<I> : E{Q) U E{Q*) —* £ be a function satisfying 

$(e*) = -l/^.(e), Ve€^(a). (17) 

Then the necessary and sufficient condition for the existence of a labeling a : ElD) C such that, 
in the notations of Fig.\^ 

$(e) = $(xo,xi) = ^ ^ $(e*) = $(yo,yi) = -^, (18) 

is given by the equations 

H ^z) = 1, H ^z*) = 1, Vxo G v{g), yo e v{g*)- (19) 

cestar(a;o;C?) c*Gstar(3/o;6*) 
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Proof. The necessity is obvious. To prove sufficiency, we construct a by assigning an arbitrary 
value (say, a = 1) to some edge of P, and then extending it successively using either of the equations 
(|18|) and the definition of labeling. Conditions (|19j) assure the consistency of this procedure. I 

Corollary 13 Integrability condition U^) for the function v : -E'(^) U E{Q*) ^ C is equivalent 
to the following one: there exists a labeling a : E{T>) C of directed edges of T>, such that, in 
notations of Fig^ 

v{yo,yi) = —( ^ = ^ • (20) 

Under this condition, the 3D consistency of the discrete Cauchy-Riemann equations is assured by 
the following values of the weights v on the edges of E{G) U E{G*) lying in the vertical faces: 

, 1 X + a 

u{y,x) = = > (21 

v\x,y) X — a 

where a = a{x,y), and X C is an arbitrary number having the interpretation of the label carried 
by all vertical edges o/D; A = a{x,x) = a{y,y). 

Proof. Apply Lemma IT^ with the function 

*(.) = (22) 

which satisfies (jlTf) due to the property Q of the weights i'. Note that in the case i^(e) G M+ the 
notation (|15|) implies that <I>(e) = exp{i(j){e)) . The formula (jlHj) with the function 1)22(1 is clearly 
equivalent to (|2U|) . To prove the second statement, we use notations of Lemma ITTl in particular 
the formula 1 + uiiiq — vim + ^oA^i = 0- According to (|2n|l . we have: vi = . Parametrize the 

(arbitrary) value of /xq as /xq = . Then it follows from the above formula that = i^^- 

An easy induction proves (|21|) for all edges in the vertical faces. I 

So, integrability of the discrete Cauchy-Riemann equations is equivalent to the existence of a 
labeling a of directed edges satisfying (|2U() . Let p : V{'D) — > C be a parallelogram realization of V 
defined by p{y) — p{x) = a(x,y). Then discrete holomorphic functions are characterized by 

f{yi) - fjyo) ^ «i - «o ^ p{yi)-p{yo) .23) 

/(a^i) - /(a^o) ai + ao p{xi)-p{xq)' 

In other words, the quotient of diagonals of the /-image of any quadrilateral {xq, yo, xi, yi) G F{V) 
is equal to the quotient of diagonals of the corresponding parallelogram. In the case of positive 
weights 1^ G R+, the labels a take values in S^, and have a geometric interpretation of edges of a 
rhombic realization of V. 

Proposition 14 The discrete Cauchy-Riemann equations \2c!\) admit a zero curvature representa- 
tion Ml\) in GL2(C)[A], with transition matrices along {x,y) € E('D) given by 

-2a{f{x) + f{y))\ 

L{y,x,a; X) = ( | , where a = p{y)—p{x). (24) 

A — a 
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Proof. This result is easy to check. It can be also systematically derived using the procedure 
outlined at the end of Sect. ID Indeed, setting A = p{x) — p{x), one writes the equation l\'2'S\i on the 
vertical face (x, y, y, x) as 

f{x) — f{y) X — a A + a f r/ \ A + a 



f{y) - fix) A + a 
where 

M{y,x,a;X) 



^ f{y) = -T-^ m + (fix) - --^ fiy)) = Miy, X, a; A) • fix), 
A — a V A — a / 



A + a (A-a)/(x)-(A + a)/(y) 
A - a 



One easily shows that these matrices form a zero curvature representation with values in GL2(C)[A], 
i.e., that holds literally, and not only projectively (up to a scalar factor). Finally, observe that 
the matrices L in (|24j) are gauge equivalent to the matrices M: 

L(y,x,a;A)= (J "^/^^) M(y, x, a; A) 

This finishes the proof. I 

The main result of the present section can be formulated as follows. Discrete Cauchy-Riemann 
equations on a quad-graph T> are integrahle if and only if they come from a parallelogram immersion 
of T> in C, weights iv being the quotients of diagonals of the corresponding parallelograms. In the 
case of real positive weights v the parallelograms are actually rhombi. 



6 Extension of discrete holomorphic functions to z 

To exploit analytic possibilities provided by 3D consistency of the discrete Cauchy-Riemann equa- 
tions, we restrict our considerations to quasicrystallic rhombic embeddings P, with the set of labels 
A = {±ai, . . . , ±ad}. Construct the two-dimensional subcomplex Jl-p in corresponding to T), as 
explained at the end of Sect. El Extend the labeling a : EiT>) — > C to all edges of Z'^, assuming that 
all edges parallel to (and directed as) carry the label a^. Now, 3D consistency of the discrete 
Cauchy-Riemann equations allows us to impose them not only on Q.u, but on the whole of Z'^. 

Definition 15 A function f : Z'^ ^ C is called discrete holomorphic, if it satisfies, on each 
elementary square of If^ , the equation 

/(n + ej +efc) - /(n) ^ a^ + a^ 
fin + Bj) - fin + efc) aj - au ' 

Obviously, for any discrete holomorphic function / : Z*^ — > C, its restriction to ViQ.x)) ~ ViJ^) 
is a discrete holomorphic function on T). To justify the reverse procedure, i.e., the extension of 
an arbitrary discrete holomorphic function on T) to Z"^, keeping the property of being discrete 
holomorphic, more thorough considerations are necessary. 

Definition 16 For a given set V C , its hull TiiV) is the minimal set Ti d IJ^ containing V 
and satisfying the condition: if three vertices of an elementary square belong to Ti, then so does the 
fourth vertex. 
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This notion is tailored for arbitrary 3D consistent four-point equations of the type Q, including 
the discrete Cauchy-Riemann equations. It is not difficult to show by induction that the hull of an 
arbitrary connected subcomplex of Z"' is a brick, i.e., a set of the type 



na,b = {n = (ni, . . . ,nrf) e Z"* : ak<nk<bk, k = l,...,d} 



(26) 



where a = (oi, . . . , a^), b = {hi, ... ,bd) are some integer vectors, with infinite values = —oo, 
bk = oo allowed. (Observe that Definition is equally well applicable to functions on bricks.) 
However, there exist combinatorial surfaces O (two-dimensional subcomplexes of Z*^), like the one 
shown on Fig. llUl that support discrete holomorphic function which cannot be extended to Z*^: the 
recursive process of extending an arbitrary discrete holomorphic function from V{^}) to its hull 
'H{V{^1)) will lead to contradictions. The reason for this is a non-monotonicity of 17: it contains 
pairs of points which cannot be connected by a path in with all edges lying in one octant. 
However, such surfaces do not come from rhombic embeddings. We will prove the absence of 
contradictions in the case of 0,x>- 



Proposition 17 For a combinatorial surface rJx" in ^ coming from a rhombic embedding of a 
quad-graph T), set 



(In case that are unbounded from below or from above on l/(J7x>)) set ak{Vtx>) = —oo, resp. 
bk{VLx>) = oo.) Then Tl^V {0,x))) = Ila,b, o.nd an arbitrary discrete holomorphic function on J^x) can 
be extended to a discrete holomorphic function on Ila,h o, unique and unambiguous way. 

For a proof of this proposition, a more detailed study of the surface f^x) will be necessary. In order 
to fix the ideas, we will assume, without loss of generality, that the circular order of the points 
±a/c on the positively oriented unit circle S-^ is the following: ai, . . . ,ad, — ai, . . . , —ad- We set 
o^k+d = ~Oik for k = 1, . . . ,d, and then define am for all m E Z by 2(i-periodicity. 

Consider the set A^, = {am, ■ ■ ■ jOm+cZ-i} of d consecutive edge slopes. The opening angle of 
the sector spanned by and am+d-i is in (0, vr). The set Am contains exactly one member ej-a^ of 
each pair ±0^, /c = 1, . . . , d. This associates to any m G the set of signs e = (ei, . . . , e^), = ±1, 
which will be denoted by e[m). The sets of signs e(m) repeat 2d-periodically, therefore not all 
possible sets of signs appear among them, but only the following 2d different ones. If m G 
then the corresponding e = e(m) is given by 




Figure 10: A non-monotone surface in 1? 



mm nfc 

neV(C-D) 



bk = bki^v) = max 

neV{C-D) 



k = 1, . . . ,d. 



(27) 




1, 1 < k < m, 
1, m < k < d. 
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and if m € [d + I, 2d], then 

+ 1, 1 < k < m - d, 
— 1, m — d < k < d. 



Fix an arbitrary xq E V{'D), and define the "sector" Um on the embedding plane C of the 
quad-graph V as the set of all points of V{'D) which can be reached from xq along paths with all 
edges from Am = {om, • • • , am+d-i}- 

This can be re-formulated in terms of f^x) as follows. Recall that the map P which identifies 
f^D with V depends on the choice of the point xq E V{'D) corresponding to E V{Qx>)- The map 
P sends Um to the set of the points of V{Qx>) which can be reached from along paths in V{^It>) 
with all edges from {eiei, . . . , e^e^}, where = efc(m) for k = 1, . . . ,d. To formulate it in a still 
another way, put into a correspondence to any set of signs e = (ei, . . . , e^) the d-dimensional octant 

Se = (eiZ+) X ... X (erfZ+) C Z*^. (28) 

In case of e = e(m), use the notation S'g(m) = Sm- Then the definition of Um is equivalent to saying 
that Um = P^^(V{il,x>) n Sm)- The following statement will be of a key importance. 

Lemma 18 The union IJm=i covers the whole of the quad- graph T>. Equivalently, the combi- 
natorial surface flx> coming from a rhombic embedding ofD lies entirely in Um=i 

Proof. Clearly, Um lies within the sector of the embedding plane with the tip at xq, spanned 
by the directions am and am+d-i- The set Am can be ordered: am -< - ■ ■ ^ am+d-i- The lower 
boundary [/~ (upper boundary C/+) of Um can be described as the path in D from the point xq 
obtained by following, at each vertex of the path, the edge with the least (resp. the largest) slope 
from Am available at this vertex, with respect to the above mentioned ordering in Am- The fact 
that V is embedded implies that all vertices of T) between C/~ and Um belong to Um- Indeed, 
suppose that there are vertices between [/~ and which cannot be reached from xq along a path 
with all edges from Am- Take such a vertex x, combinatorially nearest to xq. It cannot be reached 
from Xq along a path with the last edge from Am- Then one of the corners of one of the faces 
adjacent to x is free from edges from Am and therefore has an internal angle larger than vr, in a 
contradiction with embeddedness. Thus, Um can be described as a set of vertices between C/~ and 
f7+. Further, observe that the boundaries of the sectors Um are interlaced: Um contains all U~ 
with r E [m + 1, m + d — 2], and all U^ with r E \m — d -\- 2,m — 1]. This yields that the union of 
all C/m's covers the whole of P. ■ 
See Fig. ^2 for an illustration. 

We say that a subset T C Ila^b is an initial values locus, if, prescribing arbitrarily values of / 
on I, one can extend / in virtue of the 3D consistent Cauchy-Riemann equations in a unique and 
unambiguous way from X to the whole of Ila b (cf- | AVj ) . We will use two types of initial values 
loci. 

• Any monotone path from a to b, with all edges directed positively: 

Xi = {n.f.}^Q with no = a, njy = b, and n^+i — rij. E {ei, . . . ,ed\- (29) 

• The intersection of Ila^b with all coordinate axes: 

d 

X2 = IJ {n = nek ■ ak < n < bk}- (30) 

k=l 
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Figure 11: Sectors of the dual kagome lattice, d = 3, Ok = exp((2/c — l)7ri/6). 

Proof of Proposition 1171 We have to show that, for any set of signs e, the values of / on 
y(r2x>) n Se determine / uniquely on the hull na,b H Se- 

First, we prove this for e = e(m), so that = Sm- For the sake of notational simplicity, 
we do this for m = 1 only, i.e., for the hull TC(y{ilT>) H 5*1) = no,b- Indeed, an arbitrary point 
n G V{ilx>) n Si can be reached from along a path in V{Qx>) with all edges from {ei, . . . ,6^}. 
This is a path of the type Xi, as in (|29() . hence it is an initial value locus for the brick no,n- Since 
the union of the bricks no,n over all n € V{^lx>) n Si exhausts the brick no,b, our claim is proved. 

The bricks Ila,h ^ S^ with e ^ e(m) do not contain points of V{Qx)) in their interior. However, 
on the first step of the proof, we obtain values of / on all the coordinate axes. This gives an initial 
values locus of the type I2, as in ((SOJ, for any brick of this type. I 

Note that intersections of f^x) with bricks correspond to combinatorially convex subsets of V, as 
defined in p2] . 

7 Discrete exponential functions 

A particularly important discrete holomorphic function on Z"^ is the discrete exponential function, 
defined as 

»(-)-n(f^)'". (31) 

For d = 2, this function was considered in [FllDlj . The discrete Cauchy-Riemann equations for the 
discrete exponential function are easily checked: they are equivalent to a simple identity 

/z + aj z + ak z + ak\_aj + ak 

\z — aj z — ak J / \z — ttj z — aj — ak 

At a given n G Z'^, the discrete exponential function is rational with respect to the parameter z, 
with poles at the points eioi, . . . , edUd, where = signnfc. 
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Equivalently, one can identify the discrete exponential function by its initial values on the axes: 



e(ne,- z) = ('-±^y . (32) 

A still another characterization says that e(-; z) is the Backlund transformation of the zero solution 
of discrete Cauchy-Riemann equations on Z'^, with the "vertical" parameter z. 

Restriction of the function e(-; z) to V{Qx>) ~ ^(^) is a discrete exponential function on T> 
defined and studied in |Mll IM2[ Note that the latter depends on the choice of the point 
xq G V(T>). a question posed in ^ asks whether discrete exponential functions are dense in the 
space of discrete holomorphic functions on "D. We now show that the answer to this question is in 
affirmative, in some natural class of functions (growing not faster than exponentially). 

Theorem 19 Let f be a discrete holomorphic function on V{'D) ~ V{0,x>), satisfying 

|/(n)| <exp(C(|ni| + ... + |nd|)), yneViQv), (33) 

with some C G M. Extend it to a discrete holomorphic function on 7{{V{0,d)). Then inequality 
i'y'j^) holds for all n G 7i{V{^}Ti)), possibly with some larger constant C. There exists a function g 
defined on the disjoint union of small neighborhoods around the points ia^ G C and holomorphic 
on each one of these neighborhoods, such that 

/(n) - /(O) = ^ ^ 5(A)e(n; X)dX, Vn G nV{nv)), (34) 

where T is a collection of 2d small loops, each one running counterclockwise around one of the 
points zbofc. 

Proof. In order to extend / from V{Qx>) to TCiV^Q-z))), one makes elementary steps based on eq. 
H15|) . For instance, within the octant 5*1 these elementary steps consist of calculating the left-hand 
side of the following equation through the quantities on the right-hand side: 

/(n + e, + efe) = /(n) + ^I±^ (/(n + e,) - /(n + e,.)) • 

(In other octants everything is similar, but notations become slightly more complicated.) A simple 
induction shows that if the constant C in satisfies the inequality 



1 + 2 max 



Oj + Ok 



exp(C) < exp(2C), 



then ()33() propagates in the extension process. This proves the first statement of the theorem. 

To prove the second one, it is enough to find g{X) such that (|34|) holds on the coordinate axes, 
that is, 

/(^) - /(O) = Res 5(A) f^^y , /iti - /(O) = Res giX)(^—^Y, Vn > 0, (35) 

where fn''^ are the restrictions of / : TC{V{Qx>)) — > C to the coordinate axes: 

7^'=) = /(nefc), ak{nv)<n<bk{nv). 
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Set g{\) = Ylk=ii9k{^) + 9-k{^)), where the functions g±k{^) vanish everywhere except in small 
neighborhoods of the points ±ak, respectively, and are given there by convergent series 

«(A) = ^ (/f - /(O) + f {^y(& - &)] . (36) 

and a similar formula for 5f_fc(A). (Convergence of these series is assured by the exponential growth 
of fn^\) The easy-to-check formula 

Resi(^)" = l-(-ir, n>0, 
shows that the so defined function g satisfies (|35]). I 

8 Isomonodromic discrete logarithmic function 

We first give a construction of the discrete logarithmic function on T) which is equivalent to Kenyon's 
one jKj. This function is defined, after fixing some point xq G ViV)^ by the formula 

/(x) = ^ j|^i^e(x;A)dA, Vx G (37) 

Here the integration path T is the same as in Theorem 1191 and fixing xq is necessary for the 
definition of the discrete exponential function on T>. To make (jc{7p a valid definition, one has to 
specify which branch of log (A) is chosen around each point ±0^. This choice depends on x, and is 
done as follows. 

For each m G Z, assign to am = exp(i0m) G S"*^ a certain value of argument 9m £ choose a 
value 9i of the argument of ai arbitrarily, and then extend it according to the rule 

6m+i - Om ^ (0, tt), Vm G Z. 

Clearly, there holds 9m+d = dm + t^, and therefore also 9m+2d = + 27r. It will be convenient to 
consider the points am, supplied with the arguments 9m, as belonging to the Riemann surface A 
of the logarithmic function (a branched covering of the complex A-plane). 

The definition domain of the discrete logarithmic function is a branched covering 

CO 

U= [j Um 

m=—OD 

of the quad- graph V. Here Um is the sector Um equipped with additional data - the interval 

log{ar) e [i9m,i9m+d-i], r = m, . . . ,m + d- I (38) 

of length less than vr for the logarithms of the slopes of edges am, ■ ■ ■ , Om+d-i- If increases by 
2d, the interval on the right-hand side of (|38|) is shifted by 2Tri. Two sectors Umi and Um2 have a 
non-empty intersection, if and only if \mi — < d. It follows from Lemma ITHl that U is, indeed, a 
branched covering of V. Definition (|H7j) should be read as follows: for x G Um, the poles of e(x; A) 
are exactly the points am, ■ ■ ■ , am+d-i £ ^- Therefore, one can assume that the integration path F 
consists of d small loops around these points, and the values of log(A) at these points satisfy 
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Proposition 20 "K" The discrete logarithmic function on T>, restricted to V{Q), coincides with 
discrete Green's function on Q, up to a constant factor Itx. 

Proof. It is not difficult to see that the restriction of the discrete logarithmic function to black 
points does not branch: it is a well-defined real- valued function on V(Ci\ Clearly, this function is 
harmonic everywhere except the origin. At the origin, its Laplacian equals to the increment of / 
upon running once around the origin through its white neighbors. The values of / at the vertices 
neighboring to the origin are nothing but the arguments of the corresponding edges. Therefore, the 
above mentioned increment is equal to 27r. In order to obtain asymptotic results for the discrete 
logarithmic function, one can deform the integration path F into a connected contour lying on a 
single leaf of the Riemann surface of the logarithm, and then use standard methods of the complex 
analysis jKj. This possibility is due to the fact that functions 5^ in integral representation (|.S4|) of an 
arbitrary discrete holomorphic function, defined originally in disjoint neighborhoods of the points 
a,., in the case of the discrete logarithmic function are actually restrictions of a single analytic 
function log(A)/(2A) to these neighborhoods. I 

Now we extend the discrete logarithmic function to Z'^. To this end, recall that the sector Um of 
T) is nothing but the preimage w.r.t. P of the part of VId lying in the octant Sm C Z*^. Therefore, 
it is natural to introduce a branched covering 

00 

5* = Sya 

m=—oo 

of the set Um=i ^rn C Z'^. Here Sm is the octant Sm equipped with the set of values of log(efeafe) 
satisfying (|H5|) . Recall that = €k{m), k = 1, . . . , d, are the signs of the coordinate semi-axes of 
Sm, defined in Sect. El By definition, Smi and Sm2 intersect, if the underlying octants Sm^ and 
Sm2 have a non-empty intersection spanned by the common coordinate semi-axes, and the data 
log(efcafc) for these common semi-axes match. It is easy to see that Smi and Sm2 intersect, if and 
only if \mi — < d. 

Definition 21 The discrete logarithmic function on S is given by the formula 

/(n) = ^^i^e(n;A)dA, Vn G S, (39) 

where the integration path T consists, for n € Sm, of d loops around am, ■ ■ ■ , am+d-i; o,nd the 
branch of the logarithm on T is defined by inequality \3i^) . 

The discrete logarithmic function on T) can be described as the restriction of the discrete logarithmic 
function on 5 to a branched covering of Vtx> ~ T>. This holds for an arbitrary quasicrystallic quad- 
graph with the set of edge slopes A. 

Now we are in a position to give an alternative definition of the discrete logarithmic function. 
Clearly, it is completely characterized by its values /(nefce^) on the coordinate semi-axes of an 
arbitrary octant Sm- 

Proposition 22 For the discrete logarithmic function on S, each of d sequences fn"* = /(nefce^), 
k = 1, . . . , solves the difference equation 

n(/„+i-/„_i) = l-(-ir, (40) 

with the initial conditions 

ff^ = /(O) = 0, /f ) = fiekek) = log(efcafc). (41) 
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Explicitly, 



n 2 

/2?=E^73T' /2n+i = log(M), k = l,...,d, n>0. (42) 

£=1 

Here = efc(m), and i/ie values log(efcafc) are chosen in the interval h'j^) . 

Proof. According to eq. 1)36^ . the values fn \ with /q'^^ = 0, are defined by the expansion near 
A = efcafc, 

\ oo , 

log(A) = log(..a.) + log (-^) = /f' + Y. (^TT?) (^"« - -f"-.'- 

n=l 

This is equivalent to 

/f = logfea.-), = (44) 

The solution to these recurrent relations is given by (|42() . I 

Observe that values 1)42(1 at even (resp. odd) points imitate the behaviour of the real (resp. 
imaginary) part of the function log(A) along the semi-lines arg(A) = arg(efcafc). This can be easily 
extended to the whole of S. Restricted to black points n E 5 (those with ni + . . . + even), the 
discrete logarithmic function models the real part of the logarithm. In particular, this restricted 
function is real-valued and does not branch: its values on Sm depend on m (mod 2d) only. In 
other words, it is a well defined function on Sm- On the contrary, the discrete logarithmic function 
restricted to white points n G 5 (those with ni + . . . + re^ odd) takes purely imaginary values, and 
increases by 27rz, as m increases by 2d. Hence, this restricted function models the imaginary part 
of the logarithm. 

It turns out that recurrent relations ()4U() are characteristic for an important class of solutions of 
the discrete Cauchy-Riemann equations, namely for the isomonodromic ones. Recall the definition 
of this class. For a discrete holomorphic function / : Z*^ — >• C, the transition matrices are (cf. (EU), 

, , A + afc -2afc(/(n + efc) + /(n))\ 

Lfc(n; A) = . (45) 

\ X-Uk I 

The moving frame ^'(•, A) : GL2{C)[X] is defined by prescribing some ^(0; A), and by extend- 

ing it recurrently according to the formula 

*(n + efc;A) =Lfc(n;A)^(n;A). (46) 

Finally, define the matrices A{-; A) : Z"' ^ GL2(C)[A] by 

A(n;A) = ^5:^*-'(n;A). (47) 

These matrices are defined uniquely after fixing some A(0; A). 

Definition 23 A discrete holomorphic function / : Z*^ — > C is called isomonodromic , if, for 
some choice of A{Q; \), the matrices A{n; A) are meromorphic in X, with poles whose positions and 
orders do not depend on n G Z'^. 

^This term originates in the theory of integrable nonhnear differential equations, where it is used for solutions 
with a similar analytic characterization (1N| . 
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It is clear how to extend this definition to functions on the covering S. 
Theorem 24 The discrete logarithmic function is isomonodromic. 

This is an immediate consequence of the fohowing statement, which we formulate for functions on 
Si = for notational simplicity, but which holds actually for any octant S^. 

Proposition 25 Let 

^(0;A) = Tfn nV (48) 



A VO 0^ 

and let there he d sequences {/n*^'*}5^o satisfying, for all k = 1, . . . ,d, the recurrent relation \40j) - 
Then the discrete holomorphic function f : {l^+Y ~^ defined by the values f{nek) = fn^ on the 
coordinate semi-axes, is isomonodromic. At any point n G (^+)'^ there holds: 

A ^-^ \ \ + ai A — ai / 



=1 



with 



n ( — ^\'n,1+...+ni 

^(°)(n)= ( ^ ^ h (50) 

^0 

B<0(n) = „, I ^ -<^<"' +^^<" - 'A . C<'.(„) = n, J° + + I , (51) 

Moreover, at any point n E (Z+)'^ t/iere /loZds an isomonodromic constraint, 

d 

ni (/(n + ez) - /(n - e,)) = 1 - (-l)"i+-+-<^. (52) 

l=\ 

Proof. The matrices A satisfy a recurrent relation, which results by differentiating (|46|1 . 

A(n + e^; A) = ^^^^^L, ^(n; A) + L,(n; A)A(n; A)^-^!!; A). (53) 

Fix some k = and consider the matrices A{ne.}^\X) along the A;th coordinate semi-axis. 

Formula (|53() shows that singularities of A{ne.]^\ A) are poles at A = and at A = ia^. It is easy to 
see that the pole A = remains simple for all n > 0. As one can show (see LemmaHSlin Appendix 
IXj) . the recurrent relation (|ifl|) for /„ = /(ne^) assures that the poles A = itafc are simple for all 
n > 0. So, under condition (|4U() there holds: 

^W(nefc) , i?W(ne,) , CW(nefe) 

A(nefc; A) = h — — h — , (54) 

A A + Ofc A — afc 

i.e., eq. ()49() is valid on the /cth coordinate semi-axis, with i? (O(nefc) = CW(nefc) = for / / A:. The 
proof continues by induction, whose scheme follows filling out the hull of the coordinate semi-axes: 
each new point is of the form n + e^ ■ + e^, j 7^ A:, with three points n, n + ej and n + e^ known 
from the previous steps, where the statements of the proposition are assumed to hold. So, suppose 
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that (|49j) holds at n + e^, n + e^. The new matrix A{n + ej + e^; A) is obtained by two alternative 
formulas, 

A{n + Bj + Bk-, X) = ( ^^^^"^^ +Lk{n + Bj- X)A{n + e^; A)^ L-\n + e^; A), (55) 

and the one with interchanged roles of k and j. Eq. H55() shows that all poles of A{n + + e^; A) 
remain simple, with the possible exception of A = ±ak, whose orders might increase by 1. The 
same statement holds with k replaced by j. Therefore, all poles remain simple, and ()49() holds at 
n + ej + efc . The proof of formulas , and of the constraint is based on computations 
presented in Appendix I 

The reason for considering isomonodromic solutions on octants like is clear from ()4U() : 

indeed, this second-order ordinary difference equation has a special form, enforcing that its solution 
on the semi-axis n > is completely defined by the values at n = 0, 1, and does not depend on f_i. 
Thus, a discrete holomorphic function from Proposition 23 is uniquely defined by its initial values 
/(O) = /o and /(e,,) = /f ^ for A: = 1, . . . , d. 

Remark. The isomonodromic constraint ()52() was found in |NRGUj . without any relation 
to the discrete logarithmic function. Observe that our formulation allows us to avoid a major 
computational problem arising in |NRGUj in this context, namely that of compatibility of the 
constraint with the discrete Cauchy-Riemann equations. 

Summing up: discrete Green's function on a quasicrystallic quad-graph is the real part (i.e., 
restriction to V[Q)) of the discrete logarithmic function. The latter can be extended to a function 
on a branched covering of certain octants Sm C Z'^, m = 1, . . . , 2d. On each such octant, the discrete 
logarithmic function is discrete holomorphic, with the distinctive property of being isomonodromic. 
This function is uniquely defined either by the integral representation or by the values on the 

coordinate semi-axes \4'^ ! or else by the initial values and the constraint HSH^) . 



9 3D consistent cross-ratio equations 

The cross-ratio system is one of the simplest and at the same time one of the most fundamental 
and important nonlinear integrable systems on quad-graphs. Recall the definition of the cross-ratio 
of four complex numbers: 

. X {zq - Zl){z2 - 2:3) 

q{zo, zi,Z2, Z3) = -, (56) 

(Zi - Z2){zs - Zo) 

which yields the property 

q{zo, zi, Z2, zs) = l/qizi, Z2, ^3, 2^0)- (57) 
Let there be given a function Q : E{Q) U E{Q*) — > C satisfying the condition 

Q{t*) = 1/Q(e), Ve G E{g). (58) 



Definition 26 A function z : V{'D) C is said to solve the cross ratio equations on T> cor- 
responding to the function Q, if for any positively oriented face {xQ,yQ, xi,yi) of V there holds: 

q{z{xo), z{yo), z{xi), z{yi)) = Q{xo, xi) = l/Q(yo, yi)- (59) 
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Like in Sect. an interesting question is on the 3D consistency of the system of cross-ratio 
equations corresponding to a given function Q. 

Theorem 27 The function Q : E{g) U E{g*) C can he extended to E{G) U E{G*) giving a 3D 
consistent system of cross-ratio equations, if and only if the following condition is satisfied: 

n Q(^) n ^(^*) = 1' ^ ^(^)' ^ ^(^*)- (60) 

eestar(a;o;G) e*Gstar(yo;G*) 

Proof. We proceed as in the proof of Theorem IIUI Consider a flower of quadrilaterals around xq, 
with efc = {xQ,Xk), = {yk-i,yk)- Build the extension of this flower to the third dimension, as in 
Sect. @I Denote 

Qixo,Xk) = Q{ek) = Qk, Q{xo,yk) = l^k- (61) 

Then there holds a statement analogous to Lemma ITTI the cross-ratio equations are 3D consistent 
on the cube over the k-th. petal, if and only if 

fj.k-1 = QklJ-k- (62) 

This is donstraightforward, as in the proof of Lemma ITTI For the cube over the petal with k = 1, 
one determines on the first step the values of z at xi, yo and yi from 

q{z{xo),z{yo),z{xi),z{yi)) = Qi, 
q{z{xo),z{yo),z{yo),z{xo)) = /Uq, 
q{z{xo),ziyi),z{yi),z{xo)) = fn. 

On the second step one has three alternative values for z(xi) from 

q{z{xo),z{yo),z{xi),z{yi)) = Qi, 
q{z{yi),z{xi),z{xi),z{yi)) = ^0, 
q{z{yo),z{xi),z{xi),z{yo)) = Hi. 

A direct computation shows that these three values for z{xi) coincide if and only if /ig = Qi/^i, 
and then 

/"oyo(yi - Xq) + fj'iyiixo - yo) + Jo(yo - yi) 
fioiyi - xq) + ni{xo - yo) + {yo - yi) 

Thus, ()62() is proved. This relation yields immediately that running around xq returns back an 
(arbitrary) initial fiQ, if and only if the first condition in ()6Uj) holds. The second one follows 
similarly. I 

Corollary 28 The integrability condition I160\) for the function Q : E{g)LiE{g*) — > C is equivalent 
to the existence of a labeling c? : EiV) C of undirected edges of T>, such that, in notations of 
Fig\^ 

Q(xo,xi) = — r = (63) 

Q{yo,yi) ol{ 

This formula assures the 3D consistency of the cross-ratio equations, if one assumes that all vertical 
edges o/D carry one and the same label € C. 
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Let the labeling o? come from a labeling a : E{T>) ^ C of directed edges. Let p : 1^(1)) ^ C be 
a parallelogram realization of D defined by p{y) — p{x) = a{x,y). Then the cross-ratio equations 
are written as 

q{z{xo),z{yo),z{xi),z{yi)) = -§ = q{p{xo),p{yo),p{xi),p{yi)). (64) 

In other words, the cross-ratio of the vertices of the /-image of any quadrilateral {xo,yo, xi,yi) G 
F['D) is equal to the cross-ratio of the vertices of the corresponding parallelogram. 

Proposition 29 The cross-ratio equations \64\) admit a zero curvature representation with the 
values in GL2(C)[A], with transition matrices along {x,y) G E('D) given by 

(1 z{x) — z{y) \ 

, where a = p{y)-p{x). (65) 
AaV(z(x) - z{y)) 1 j 

Proof. This result is easy to check. To derive it systematically using 3D consistency and the 
procedure outlined at the end of Sect. |1J observe that eq. H64() on the vertical face (x, y, y, x) can 
be written as 

q{z{x),z{y),z{y),z{x)) = W 4^ z{y) - z{y) = L{y,x,a; X) ■ {z{x) - z{x)), 

where A = {p{x) — p(x))~^, and the matrices L{y, x, a; A) are as in One easily shows that these 
matrices L form a zero curvature representation with values in GL2(C)[A], i.e., that (fTT|) holds as 
it stands, and not only up to a scalar factor. I 

The main result of the present section can be formulated as follows. Integrable cross-ratio 
equations on a quad-graph T> come from parallelogram immersions of D in C, the coefficients Q 
being the cross-ratios of the corresponding parallelograms. In the case of unitary values Q £ the 
parallelograms are actually rhombi. 



10 Circle patterns and the cross-ratio system 

Definition 30 A Delaunay decomposition of C is a cell decomposition Q such that the boundary 
of each face is a polygon inscribed in a circle, and these circles have no vertices in their interior. 
These circles build a circle pattern with the combinatorics of Q. 

The vertices z : V{Q) ^ C of a Delaunay decomposition are the intersection points of the circles 
of the corresponding pattern. The circle of the pattern corresponding to a face y E F{Q) will be 
denoted by C{y). If two faces yo,yi G F{Q) have a common edge (xo,xi), then the circles C(yo) 
and C{yi) intersect in the points z{xq),z{xi). In other words, the edges of Q correspond to pairs 
of neighboring (intersecting) circles of the pattern. Similarly, if several faces yi,y2, ■ ■ ■ ,ym of G 
meet at one vertex xq G V{Q), then the corresponding circles C(yi), C{y2), . . . , C{y„i) also have a 
common intersection point z{xq). 

Given a circle pattern with the combinatorics of ^, we can extend the function z to the vertices 
of the dual graph Q* , setting 

z{y) = center of the circle C(y), y G F{Q) ~ V{Q*). 

After this extension, the map z is defined on all of V{T>) = V{Q) U V{Q*), where T> is the double of 
Q. Consider a face of the double. Its vertices xq, yoi Ui correspond to the intersection points and 
to the centers of two neighboring circles Co, Ci of the pattern. The following statement is obtained 
by a simple computation. 
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Figure 12: Circle pattern 



Lemma 31 If (j) is the intersection angle ofCo,Ci, as on Fia M!A then 

q{z{xo), z{yo), z{xi), z{yi)) = exp{2i^). (66) 




Figure 13: Two intersecting circles 



It will be convenient to assign the intersection angle cp of C{yo), C{yi) to the edge {yo,yi) € E{Q*). 
Extend the function : E{g*) (0,7r) to E{g) by setting 0(e) = vr - (^(e*). 

Proposition 32 Let Q he Delaunay decomposition of a plane, and consider a circle pattern with 
the combinatorics of Q and with the intersection angles (p : E{Q*) — > (0,7r). Let {z{x) : x G V^Q)} 
and {z[y) : y € V{Q*)} consist of intersection points of the circles, resp. of their centers. Then 
z : V{T>) —>■ C satisfies a system of cross-ratio equations with the function Q : E{Q) U E{Q*) §^ 
defined as Q{t) = exp(2i0(e)). There holds: 

J] exp(2i(/)(e)) = 1, Vxo G ViG). (67) 

eGstar(a;o;6) 
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The following condition is necessary and sufficient for the integrability of the system of cross-ratio 
equations: 

H exp(2i<A(e*)) = 1, VyoG^r), (68) 

c*Gstar(j/o;6*) 

i.e., for each circle of the pattern the sum of its intersection angles with all neighboring circles of 
the pattern vanishes (mod vr). 

Proof. The relation (|67|) is obvious for geometrical reasons: for an arbitrary common intersection 
point of circles of the pattern, the sum of their consecutive pairwise intersection angles vanishes 
(mod vr). Now the claim follows from Theorem 1271 I 

We can formulate the main result of this section as follows. Combinatorial data Q and inter- 
section angles (p ■ E[Q*) (0, vr) belong to an integrable circle pattern, if and only if they admit 
an isoradial realization. This latter realization gives a rhombic immersion of the double T), and 
generates also a dual isoradial circle pattern with the combinatorial data Q* and intersection angles 

11 Hirota system 

We have seen that integrable circle patterns deliver solutions z : V{T>) — > C to integrable cross- 
ratio systems. These solutions are characterized as follows: the z-image of any quadrilateral 
ixQ,yo,xi,yi) from F{T>) is a kite with two pairs of sides of equal length (incident to the white 
vertices -z(yo)) -^^(yi))) with the prescribed angle vr — (/> at the black vertices z{xq), z{xi). The 
following transformation of the cross-ratio system is useful in order to single out this class of kite 
solutions. 

Definition 33 Let a : E{T>) — > C 6e a labeling, and let p : V{T>) ^ C be the corresponding 
parallelogram realization ofV defined by p{y) —p{x) = a{x,y). A function w : V{T)) —>■ C is said to 
solve the corresponding Hirota system, if for any positively oriented face {xo,yQ,xi,yi) € F{T>) 
there holds, in the notations of Fig^ 

aow{xo)w{yo) + aiw{yo)w{xi) - aow{xi)w{yi) - aiw{yi)w{xo) = 0, (69) 

or, in a more invariant fashion, 

w{xo)w{yo){p{yo) - p{xo)) w{yo)w{xi){p{xi) -p(yo)) + 

w{xi)w{yi){p{yi) - p{xi)) + w{yi)w{xo) {p{xo) - p{yi)) = 0. (70) 

Obviously, a black-white scaling, i.e., a transformation w ^ cw on V{Q) and w c~^w on V{Q*) 
with a constant c, maps solutions of the Hirota system into solutions. We will identify solutions 
related by such a transformation. 

Proposition 34 Let w : V{T>) ^ C be a solution of the Hirota system. Then the relation 

z{y) - z{x) = a{x,y)w{x)w{y) =w{x)w{y){p{y) -p{x)), \f{x,y) e E{T>), (71) 

correctly defines a unique (up to an additive constant) function z : V{T>) — > C which is a solution of 
the cross-ratio system ( |6'^| ). Conversely, for any solution z of the cross-ratio system ( |6'^| ), relation 
l |yi| ) defines a function w correctly and uniquely (up to a black-white scaling); this function w solves 
the Hirota system Hf)fJ\). 
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Proof. Simple calculation based on closing conditions around the quadrilateral {xQ,yQ,xi,yi). I 



Proposition 35 Let all a € S^, so that p : V{'D) C is a rhombic realization ofT>. Let z : 
V{'D) C be a solution of the corresponding cross-ratio system \64\) - It corresponds to a circle 
pattern, if and only if the corresponding function w satisfies the condition 

w{x)eS^, w{y)en+, ^xeV{g), yeV{g*). (72) 

The values w : V{Q*) ]R_|- have then the interpretation of the radii of the circles. 

Proof. The function z corresponds to a circle pattern, if and only if all elementary quadrilaterals 
{z{xo), z{yo), z{xi), z{yi)) are of the kite form with the properties: 

- the pairs of edges incident with white vertices have equal length, 

- the angles at black vertices are equal to the corresponding angles of the underlying rhombi. 
As easily seen, these conditions are equivalent to: 

\w{xo)\_^ and ^GM„ 



\w{xi)\ w{yi) 

respectively. This yields (|7^ . possibly upon a black- white scaling. I 

Remark. The conditions H72|) form an admissible reduction of the Hirota system corresponding 
to a rhombic realization, in the following sense: if any three of the four points w{xq), w{yQ), w{xi), 
w{yi) satisfy the condition ()72() . then so does the fourth one. This is immediately seen, if one 
rewrites the Hirota equation (|69|1 in one of the equivalent forms: 

w{xi) _ aiwjyi) - aoUj{yo) ^ w{yi) _ apw^xo) + aiw{xi) ^^^^ 
w{xo) aiw{yo) - aow{yi) w{yo) aow{xi) + aiw{xo)' 

Proposition 36 a) Let a : E{T>) <C he a labeling, and let p : V{T>) ^ C be the corresponding 
parallelogram realization of D defined by p{y) — p{x) = a{x,y). Then the corresponding Hirota 
system is 3D consistent. 

b) Let all a G S^ so that p : C is a rhombic realization of T>. Consider a solution 

w : V{T>) — > C corresponding to a circle pattern with the combinatorics of Q, i.e., satisfying \7i^ . 
Consider its Bdcklund transformation w : V{'D) C with an arbitrary parameter A G §^ and with 
an arbitrary initial value w{xo) G M+ or id{yo) G S^. Then there holds: 

wix) G R+, iviy) G §\ Vx G ViG), y G ViG*), (74) 

so that w corresponds to a circle pattern with the combinatorics of G* • 

Proof. Statement a) is a matter of a direct computation. In the notations of Fig.lHl suppose 
that the Hirota equation ()69() holds on all faces of the cube, wherein the vertical edges carry the 
(arbitrary) label A. One finds that all three alternative ways to compute vuixi) lead to one and the 
same result, namely 



_ A(ag - al)w{yo)w{yi) + ai(A^ - al)w{yo)w{xo) + aojal - X'^)w{y i)w{xo) 
X{al - af )w(xo) + ai(A2 - ag)w(yi) + Qo(af - X^)w{yo) 

Statement b) follows from the Remark above. 
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Proposition 37 The Hirota system h6y\) admits a zero curvature representation with the values in 
GL2(C)[A], with transition matrices along {x,y) € E{'D) given by 

1 —aw{y) 

L{y,x,a;X)= | | , where a = p{y)—p{x). (75) 

-Xa/w{x) w{y)/w{x) 

Proof. The matrix (|75)) is obtained also directly from (|65)) by the substitution (|71jl . followed by a 
simple gauge transformation 

(I \ /I 

Alternatively, a systematic derivation of this result is based on the 3D consistency and the procedure 
outlined at the end of Sect. |3 I 

The main result of the present section is as follows. Integrahle circle patterns can he alternatively 
described by solutions of the Hirota system with a special property of being real-valued on V{Q) and 
unimodular on V{Q*). 

12 Linearization 

Let Q : ^ C be a labeling, and let p : V(T>) — > C be the corresponding parallelogram 

realization of P defined by p{y) — p{x) = a{x, y). Then the formula 

zo{x) = p{x), wo{x) = 1, Vx G V{V), 

gives a (trivial) solution of the cross-ratio system ()64() and the corresponding (trivial) solution of 
the Hirota system. Suppose that zo : V{'D) — > C belongs to a differentiable one-parameter family of 
solutions Ze : V{'D) — > C, e € (— eo, eo), of the same cross-ratio system, and denote by w^ : V{'D) — > C 
the corresponding solutions of the Hirota system. Denote 



dz^ 



6=0 V de y 



Proposition 38 Both functions /, g : V{T>) C solve discrete Cauchy-Riemann equations h2'J\) . 
Proof. By differentiating H71j) . we obtain a relation between the functions f,g : V{T)) — > C: 

g{y)-9{x) = {f{x)+f{y)){p{y)-p{x)), y{x,y) e E{V). (77) 

The proof of proposition is based on this relation solely. Indeed, the closeness condition for the 
form on the right-hand side reads: 

(/(xo) + f{yo)) {p{yo) - p{xo)) + (/(yo) + f{xi)) {p{xi) - p{yo)) + 

{f{xi) + f{yi)) {p{yi) - p{xi)) + (/(yi) + /(xq)) (p(xo) - p{yi)) = 0, 



which is equivalent to (|23j) for the function /. Similarly, the closeness condition for /, that is, 
{f{xo) + f{yo)) - (/(yo) + /(xi)) + (/(xi) + fiyi)) - {f{yi) + /(xo)) = 0, 
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yields: 



9{yo) -g{xo) _ g{xi) - gjyo) ^ g{yi) - gjxi) _ gjxo) - g{yi) 

p{yo)-p{xo) p{xi)-v{yo) p{yi)-p{xi) p{xo)-p{yi) 



= 0. 



Under the condition p(yo) — p{xq) = p{xi) — p{yi), this is equivalent to (|23() for g. I 

Remark. This proof shows that, given a discrete holomorphic function / : 1^(1') — C, relation 
H77() correctly defines a unique, up to an additive constant, function g : V{'D) — > C, which is also 
discrete holomorphic. Conversely, for any g satisfying the discrete Cauchy-Riemann equations ()23() . 
relation (|77|) defines a function / correctly and uniquely (up to an additive black- white constant); 
this function / also solves the discrete Cauchy-Riemann equations (j^SJ- Actually, formula ((77|) 
expresses that the discrete holomorphic function / is the discrete derivative of g, so that g is 
obtained from / by discrete integration. This operation was considered in |DH ID2| IMlj . 

Summarizing, we have the following statement. 

Theorem Z9 a) A tangent space to the set of solutions of an integrable cross-ratio system, at 
a point corresponding to a rhombic embedding of a quad-graph, consists of discrete holomorphic 
functions on this embedding. This holds in both descriptions of the above set: in terms of variables 
z satisfying the cross-ratio equations, and in terms of variables w satisfying the Hirota equations. 
The corresponding two descriptions of the tangent space are related by taking the discrete derivative 
(resp. anti- derivative) of discrete holomorphic functions. 

b) A tangent space to the set of integrable circle patterns of a given combinatorics, at a point 
corresponding to an isoradial pattern, consists of discrete holomorphic functions on the rhombic 
embedding of the corresponding quad- graph, which take real values at white vertices and purely 
imaginary values at black ones. This holds in the description of circle patterns in terms of circle 
radii and rotation angles at intersection points (Hirota equations). 

13 Isomonodromic discrete power function 

Like in Sect. |21 one can consider functions z : Z*^ — > C and : Z'* — > C, satisfying, on each 
elementary square of U^, the cross-ratio and the Hirota equation, respectively, and ask about 
isomonodromic solutions. As shown in |ABll IAB21 IBHj . this leads to a discrete analog of the 
power function. Since the latter references contain a detailed presentation of these results in terms 
of the cross-ratio variables z, we restrict ourselves here to similar results in terms of the Hirota 
variables w. (Recall that transition matrices in these two formulations actually coincide, up to a 
simple gauge transformation by diagonal matrices which do not depend on A.) 
Transition matrices for the Hirota system on are: 



With these transition matrices, isomonodromic solutions are defined in exactly the same manner 
as in Sect. El 

Proposition 40 Let 



and let there be d sequences {wn }^o satisfying, for all k = 1, . . . ,d, the recurrent relation 




(78) 




(79) 




(80) 
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Then the solution w : (2+)'^ C of the Hirota system, defined by the values w{nek) = Wn^ on the 
coordinate semi-axes, is isomonodromic. At any point n G there holds: 

with 

AW(n)= I "; 1 , (82) 




w{n + ei) + w{n-ei) 1 i/q,^ w(n - e,) 



(83) 

The upper right entry of the matrix A^^\yi), denoted by the asterisk in is actually given by 



4?(n) = -E^SH- (84) 



Z=l 

Moreover, at any point n G (^+)"' t/iere /loZds an isomonodromic constraint, 
d 

in J 

w{n + ei) +w{n- ei^ 

Proof. The scheme of the proof is the same as for Proposition 1251 Fix some k = l,...,d, 
and consider the matrices j4(nefc;A) along the A;th coordinate semi-axis. It follows from formula 
H53() that singularities of A{nek; X) are poles at A = and at A = a^^- While the pole A = 
automatically remains simple for all n > 0, this is not necessarily so for the pole A = a^^- As one 
can show (see Lemma IIH in Appendix^l); the recurrent relation ()8U() for Wn = f{nek) assures that 
the pole A = is simple for all n > 0, and 

A{nek; A) = ^ + ^ ^ , (86) 

A A -a;. 

i.e., eq. H81() is valid on the kth. coordinate semi-axis, with B^^\nek) = for I ^ k. To prove that 
eq. (|8T|) is valid also elsewhere, one argues by induction: suppose that (|8T|) holds at n-|-ej, n-|-efc. 
Then eq. H55() shows that all poles of j4(n + ej -|- e^; A) remain simple, with the possible exception 
of A = , whose order might increase by 1. The same statement holds with k replaced by j. 
Hence, all poles remain simple. Therefore, H81() holds at n + + e^, possibly up to a term which 
does not vanish with A oo. Such a term is absent, if the right-hand side of (|53() vanishes with 
A oo, that is, if 

])(AW(n)+X:i?(')(n))(° [J)=0. (87) 

Clearly, the latter equation is equivalent to (jHlJ . Computations towards the proof of ()84|) , as well 
as of (|82|) . (|83j) and of constraint (|85|) . are presented in Appendix IbI I 

Remark. Again, the isomonodromic constraint ()85() was found in |NRGUj . In the approach of 
that paper, consistency of the constraint with the Hirota equations (called lattice MKdV there) is a 
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difficult problem, only manageable with the help of a computer system for symbolic manipulations. 
In our formulation, this comes for free, as a natural consequence of the construction based on the 
3D consistency of the Hirota system. 

A solution of the Hirota system given in PropositionUHlis completely defined by its initial values 
w{0) = Wo and w{ek) = w^^^ iov k = 1, . . . ,d. The choice 

Wo = 1, ^ = exp(ipfc) , k = l,...,d, 

with arbitrary constants pk, leads to the following solution on the semi-axes: 



2n=n^7^' -?ii=exp(^P.). (89) 



Observe the asymptotics at n ^ oo, 



w 



(fc) 

2n 



3(7)n2T-i(l + 0(n-i)). (90) 



The following special choice of pk defines the discrete analog of the function w i-^ w'^'^ ^ on (Z+)'^: 

ipA: = (27 - l)logaA,., so that w{ek) = al'^~^ . (91) 

Remark. In the variables z the initial values {zii^}^^^ on the semi-axes are given by the 
following analog (and consequence) of H80|): 

^ i^n+l-Zn)iZn-Zn-l) ^ ^^^^ ^g^) 
Zn+1 Zji—l 

and then the corresponding solution of the cross-ratio equation satisfies an analog of (|85|) : 



d 



(^^(n + ej) - z(n))(z(n) - zjn - ej)) 

/ '^j 7 ^ } ^ = iz{n). (93) 

fr{ z{n + ej)-z{n-ej) 

Recall also |ABH IXB 2 BH that the discrete analog of the function z z^'^ on (Z_|_)'^ is charac- 
terized by the constraint ()92() and the following choice of the initial conditions: 

z(0)=0, z(efc) = a27^ k = \,...,d. (94) 

Clearly, this choice of initial conditions is equivalent to (|88j) . if one takes into account the basic 
relation 1)711) between the variables z and w. 

Like in Sect. |SJ isomonodromic solutions similar to those of Proposition ^U] can be defined not 
only on (Z+)'^ but on any octant S^. They are characterized by the initial data 

t/;(0)=0, Z£;(efcefc) = (efcQfc)2^-\ A; = l,...,d, (95) 

and give discrete analogs of the function w ^ w'^'^^^ on ^g. Such a solution is fixed by an inde- 
pendent choice of branches of the function w'^^~^ at the points w = e^a^. This is equivalent to 
choosing the branches of the function logw, because of w'^^~^ = exp((27 — 1) logw). 

Definition 41 The discrete power function w'^'^^^ on S is a complex-valued function whose 
restriction to Sm is defined as the unique isomonodromic solution w : Sm — > C 0/ the Hirota system 
on the corresponding Sm with the initial data i9,5)) fixed by the condition 
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Clearly, the discrete power function takes real values at the white points and unimodular values at 
the black points, so that it corresponds to a circle pattern. 

Proposition 42 The tangent vector to the space of integrable circle patterns along the curve con- 
sisting of patterns w'^^~^ , at the point corresponding to j = 1/2, is the discrete logarithmic function. 

Proof. We have to prove that the discrete logarithm / and the discrete power function w'^'^~^ are 
related by 

It is enough to prove this on the coordinate semi-axes of each octant Sm- But this follows by 
differentiating with respect to 7 constraint ()8U() and initial conditions (jUSJ) at the point 7 = 1/2, 
where all w = 1: the results coincide with (|4()|) and (|41j) . respectively. I 



14 Concluding remarks 

Results of Sect. ^Jcan be generalized to the case of linearization at an arbitrary (not necessarily 
parallelogram) solution zq of the cross-ratio system and the corresponding solution vjq of the Hirota 
system. In this case the relation between derivatives ((7H|) . coming to replace eq. ((77j) . reads: 

9{y)-9{x) = {f{x)+f{y)){zo{y)-zo{x)), y{x,y) e E{V). (96) 

Arguments similar to those of the proof of Proposition |2H1 show that in this case the function / is 
discrete holomorphic with respect to zq, i.e., 

fjyi) - fjyo) zo{yi) - zojyo) 

-jr~\ — 7r^ = ~f — ^ r^' y{xo,yo,xi,yi) e f{v). 

J[xi)-f{xo) zo{xi) - zo{xo) 

The function g, in general, is no longer discrete holomorphic. 

Thus, a tangent space to the set of integrable circle patterns of a given combinatorics, at an 
arbitrary point, consists of functions, discrete holomorphic with respect to the kite-form embedding 
Zq of the corresponding quad-graph. This holds for the description of circle patterns in terms of circle 
radii and rotation angles at intersection points (Hirota equations). The elements of the tangent 
space are characterized by the property of being real at white vertices and purely imaginary at 
black ones. 

A number of constructions of the present paper can be generalized to the case of kite-form 
(rather than rhombic) embeddings coming from an integrable circle pattern. In particular, differ- 
entiating the discrete vlP''^~^ with respect to 7 at a point 7 7^ 1/2, one obtains a sort of the discrete 
logarithmic (and Green's) functions on the kite- form quad-graph corresponding to z^'^ . 



A Appendix: proof of Proposition 1251 

Lemma 43 Let the matrix A(0; A) he as in J^<$| ). Fix some k = 1, . . . ,d. Then singularities of the 
matrices A{nei.; X) are poles at X = 0, X = ia^. For n > 0, the poles A = and X = —a^ are 
simple. The pole X = Ok is simple for all n > 0, if and only if recurrent relation i40[ ) holds for 
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fn = f{nek). In this case there holds \54\l with 

yl(0)(ne,)= T ' 1 , (97) 






Proof proceeds by induction. Putting matrices (|45j) into the recurrent definition (|5.S|) . one finds 
immediately: 

Tl Ti 

Aii{nek] A) = — , ^22(nefc; A) 



\ + ak A - afc 

and the following recurrent relation for the upper right entry of the matrix A{-; A), 

^l2((^^ + Ijefc; A) = ^12(^0^; A) + ' 



A - Ofc A - afc V A + Ofe A - Ofc 

Assume that there holds eq. ()54|) . Then it holds also with n 1-^ n + 1, if and only if no higher order 
pole appears at A = afc by this transition, what is equivalent to 

&S{nek)=n{fn+i + fn). (99) 

This has to be considered as a recursive definition of the sequence {fn} (isomonodromic constraint). 
Under this condition, we find the recurrent relations for the upper right entries of the matrices A^^^ , 

4°^((n + l)efc) = -4?(nefc), (100) 
5g)((n + l)e,.) = -(n + !)(/„+!+/„), (101) 
Cg)((n + l)efc) = (n + l)(/„+i + /„)+i?i^)(nefc) + cg)(nefc) + 24?(nefc). (102) 

Now (Unni, Itrnri) yleld: 

Af^inek) = i-ir, B[l\nek) = -n(/„ + U-i). (103) 
Adding all three equations ()100|) - H102() . we find that 

a'S {nek) + i?2^ (nek) + ci^^ (ne,.) = 1, (104) 

and this together with (|99|) . (|103jl implies explicit form (|40p of the isomonodromic constraint. I 

Proof of Proposition 1251 continued. As shown in the main text, the induction from n + e^, 
n + e/j to n + + proves formula (|49|) . From the diagonal part of eq. (|55|) one easily derives 
that for all n G Z'^, 

d d 

ni A ^ sr^ '^l 



An(n;A) = ^, ^22(11; A) = ^ 



A + a/ ^-^ A — a; 

1=1 ' /=i ' 

The following formula is an easy consequence of eq. under the limit A ^ 00: 

d 



1=1 
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It remains to show that the following relations propagate in the evolution defined by the recurrent 
relation (f^ : 

45(n) = (106) 
Bg(n) = -nK/(n) + /(n - e,)), cg(n) = n,(/(n + e,) + /(n)). (107) 

Indeed, constraint follows then immediately, because it coincides with H105|l . if one takes (|106() - 
(jlOTjl into account. Writing now the upper right entry of eq. in length, we find the following 
recurrent relations: 

4°2(n + ej+efc) = -4?(n + ej), (108) 

^12 (n + + Sfc) = , ^12 (n + ej) ^7^(?^/ + Sij + 6ik){f{n + + e^) + /(n + e^)), 

a« + Ok ai + Qfc 

(109) 

C7g(n + e,+efc) = ^i±^ C7g(n + e,-) - ^^(n; + <5,,-)(/(n + e, + e,) + /(n + e,)), 

a« — afc a/ - «fc 

(110) 

the latter formula being valid for / ^ k only. For c['^{n + ej + e^) there holds a similar but 
much longer formula, which we actually will not need. Now, eq. H108|) readily yields H106|) . By the 
induction hypothesis, eqs. ()109() and (|110() with / ^ k can be rewritten as 

B[2in + ej + efc) = 

-{ni + 6ij) (^i^{f{n + e,) + /(n + e, - e,)) + ^^(/(n + + e^) + /(n + e,))) , 
Cf^{n + ej + ek) = 



{ni + 6ij) 



(/(„ + ej + eO + /(n + e,)) - ^^(/(n + e, + e^) + /(n + e,))) . 



But the discrete Cauchy-Riemann equation for the corresponding elementary squares imply that 
the latter two equations are equivalent to 

B[^^{n + ej +ek) = -{ni + 6ij){f{n + ej + ek) + f{n + ej + ek - ei)), 
C[2{n + ej + ek) = {rii + 6ij){f{n + ej + + ei) + f{n + ej + ej,)), 

which coincide with (jlOTj) at n + ej + for I ^ k. By interchanging the roles of k and j, we see 
that p()7|) at n + ej + holds also for / / j, and thus for all / = 1, . . . , d. I 



B Appendix: proof of Proposition 1401 

Lemma 44 Let the matrix A{0; A) be as in ^7^. Fix some k = 1, . . . ,d. Then singularities of the 
matrices A{nek; A) are poles at X = 0, X = . The pole X = is simple. The pole X = ol^ is 
simple for all n > 0, if recurrent relation Hy^) holds for Wn = w{nek). In this case there holds \8(^) 
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with 

A(^\nek)= I 1 , (111) 





SW(nefc) = ■ . (112) 

Wn+l+Wn-l yi/ak Wn-l J 

The upper right entry of the matrix A^'^^nek), denoted in by the asterisk, is given by 

Af^{nek) = -B[l\nek). (113) 

Proof. Assume that eq. (|86() holds. Put matrices H78|) into recurrent definition (|53|) . It is easy to 
see that no higher order pole appears at A = by the transition n i-^ n + 1, if and only if 

1 -'"^''-A Bi^Hne,) ( ' 1=0. 

-l/{akWn) Wn+l/WnJ \ 1/ (Ofclfn+l ) Wn/Wn+1 ^ 

This is equivalent to 

(l -akWn+i) B^'^^nek) ( ^ J = 0, (114) 

or, written in length, to 

Bf^{nek) - BlWnek) + akWn+iB^Wnek) B['^\nek) = 0. (115) 

This is a recursive definition of the sequence {wn} (an isomonodromic constraint). Notice that this 
is a quadratic equation for Wn+i, unlike (|99() . which was a linear equation for fn+i- 

Under condition H114|) . or ()115() . eq. H86|) holds also by n i— > n + 1, possibly with an additional 
A-independent term on the right-hand side, which vanishes if and only if ()113|) holds. We will show 
in a moment that this is indeed the case. One readily finds recurrent relations for the matrices 

(nek) and B'^^\nek). For the matrix A(°)(nefc) they read: 

.rnw/ 1N X -akWn+i\ ,/l akWn \ 

A^'^>{{n + l)ek)= \ \A^'^\nek)\ \. (116) 

Wn+l/Wnj Wn/Wn+lj 

This proves formula Ulllj) . with a recurrent relation for the upper right entry: 

4?((n + l)efc) = -^akWn + ^4?(nefc). (117) 

Wn+l 

For the matrix B^^\nek) the recurrent relations read, in components: 

5i?((n + l)efc) = -f + Bg\nek) ^ (4°2 (ne^) + ^{^(nefc)) , (118) 

Bg)((n + l)efc) = l-^ + i3j^)(nefe) + ^ {A^^ (nek) + B['^\nek)) , (119) 

Bg)((n + l)efc) = ^akWn-^^Af^inek), (120) 

Bg)((n + l)efc) = ^ + ^ (4?(nefc) + 2i?i^)(nefc)) 

^ 'Bi'^{nek)-B['l\nek)). (121) 
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Comparing now H12U|) with (jllTj) . we see that H113|) holds for all n, as claimed above. Upon using 
this fact and constraint (|115j) . we can rewrite formulas (|11<S|) - (|12H) as follows: 

B['l\{n + l)ek) = ^ + Bg\nek), (122) 
Bg\{n + l)ek) = l-^z + Bifinek), (123) 
5g^((n + l)efc) = ^akWn + ^^B[l\nek), (124) 

Wn+l 

<((n + l)e,) = i^ + !^i?g)(ne,). (125) 

akWn Wn 

These relations together with constraint (|115() define the sequence {wn} and the matrices B^^\nek) 
completely. First of all, there follows from (|122j) . (|123|) : 

Bf^{nek) + B^^{nek) = n, B^^ {n^k) - B^^ {nek) = (t " (l " (-1)")- (126) 
Further, there follows from (fT^ - (fn5| ): 

B^^{{n + l)ek)-^—Bf^{{n + l)ek) = Bg\nek) ^('^(nefe), (127) 

akWn akWn+1 

{{n + l)ek) - akWnB^^i ii^ + '^)^k) = {nek) - akWn+iB'^i {nek). (128) 

Subtracting these two equations and taking (|115j) into account, we find, upon downshifting n: 

Bf^{nek) - Bg\nek) + afcu;„_iS^f (ne^) B[l\nek) = 0. (129) 

This yields that one of the solutions of eq. (jllSp . considered as a quadratic equation for Wn+i, is 
Wn+i = —Wn-1- We will be interested in the second one. To find it, add eqs. (|115|) and (|129|) and 
derive, under the condition Wn+i + Wn-i 7^ 0: 

B[2'{nek) = alwn+iWn-iB^Wnek). (130) 

Due to H115() . the right-hand sides of (|127() . H128() are equal to one another. Using there (|130() . we 
finally come to 

Bu{nek) = akWn+iB!2'l\nek), B^l\nek) = akWn^iB^;^^ {nek). (131) 
This together with (|126|) yields 

B^hn^k) = —, r, (132) 

ak{Wn+l + Wn-l) 

and now both the expression pi2|) and explicit form (|8()j) of the constraint (|115j) follow readily. I 

Proof of Proposition I40L continued. As demonstrated in the main text, if eq. ()81|) holds 
at n + Gj, n + e^, then it holds also at n + ej + e^, provided eq. (|51)) is valid. To prove eq. (|51)) . 
put expression (|78() into eq. H55() . and find recurrent relations for the matrices ^('^^(n) and B^^\-n). 
Upon use of the abbreviation Lj^k{^) = Lk{n + ej; A), we have: 

A(0)(n + e, +efe) = L,- ^(O)^^") (n + e,)L,+fc,_fc(0), (133) 

S(')(n + ej- +efc) = Lj^k{ar^) , o J -^j+fc.-fc(«r^)' ^ ^ ^' 

1 - afoi 

B('\n + e,+ek) = -L,,,(a,-^) f ^^(n + e,) + ^ ^^'^^ V ) ^.+fe,-feK'') 

+ lower triangular matrix. (135) 



(134) 
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Here we used the fact that L~^(A) = Lj+fc^_fc(A)/(l — Aa^^). Taking into account that the upper 
triangular part of the matrix Lj^k{\) does not depend on A (and coincides with Lj^fc(O)), we see 
that eqs. (|133j) - (|135j) imply the desired property (|8l|) . which proves formula (|8T|) at n + + e^. 
After eq. (j5T|) has been proved, it is instructive to rewrite eq. as 

^(n + e^; A)Lfc(n; A) - Lfc(n; A)^(n; A) = 
and consider the limit A — > cxo of this formula. Due to ()78() and ()81|). this limit reads: 

^(°)(n + e.) + X:i^(0(„ + e,.))(; - {^^ q) (^^^n) + X; i?«(n)) = °) . (136) 

Clearly, this equation contains more information than 1)8 7|) . More precisely, the diagonal terms of 
(|136|) are equivalent to 1)8 7|) . while the lower left entry gives an additional identity. 

Further, eq. (|133j) yields immediately that the matrix A^^\n. + + e^) retains the upper 
triangular form 1)821) . Turning to the matrices B^^\yy), observe first of all that formula H83|) is 
equivalent to 



i?«(n) ~ (l a,«;(n-e;)) , tri?«(„) =nz, (137) 




where the sign ~ means "equal up to a scalar factor". Suppose that this holds at the points 
n + Gj and n + for all / = 1, . . . , d. Then it follows from ()134|) that for all / 7^ A; the matrix 
^{0(n + ej + e/j) is also of rank 1, and its trace is equal to tr i?(')(n + e^) = n/ + 5ji = ni + 6ji + 6m- 
It remains to prove that 

(1 a;?i;(n + ej - e;)) Lj+fc -fc(a;"^) ~ (1 a/u;(n + + - e/)) . (139) 
These equations, written in length, read: 

wn + ej ^ — ^ , — ■ — r = w[n + ej + ek + ei), (140) 

akw{n + ej + e,) - aiw{n + ej + e^) 

w{n + ej) — — — ■ — = w(n + ej + ek-ei), (141) 

akw{n + Bj - ei) + aiw{n + ej + Bk) 

and are nothing but the Hirota equations on the corresponding elementary squares. Thus, recurrent 
relation H134|) implies that formula ()137() holds at n + ej + for I ^ k. By interchanging the roles 
of j and k, formula (|137|) holds for I ^ j, and thus for all / = 1, . . . , d. It remains to prove the 
isomonodromic constraint (|85|) . But it is not difficult to see that it is a direct consequence of the 
lower left entry of eq. (|136|) , if one takes into account expressions (|82|) and (|83|1 . I 
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